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THE MODEL 

The Nambu-Jona-Lasinio Model extended to include the 
't Hooft determinant in the light quark sector (u,d and 
s) is regarded as a useful tool for the study of low en- 
ergy hadron phenomenology as it shares with QCD the 
global symmetries, incorporates by construction a mech- 
anism for dynamical chiral symmetry breaking and in- 
cludes explicit breaking of the Ua(1) symmetry. It has 
been shown that the inclusion of eight quark interactions 
can solve a fundamental problem of this model: the ab- 
sence of a globally stable ground state [1]. Several new 
features emerge from the addition of the 8q terms and 
are to a large extent dictated by the part of these interac- 
tions which violates the Okubo-Zweig-Iizuka (OZI) rule. 
At the same time they leave open the possibility of al- 
most reproducing the same mesonic stectra [2] that are 
obtained without them, thus preserving the validity of 
the model. It was shown that the critical endpoint in the 
phase diagram of this model is pushed towards vanishing 
chemical potential and higher temperatures with increas- 
ing strength of the OZI-violating eight-quark interactions 
[3]. The introduction of the eight quark interactions also 
affects the number of effective quark degrees of freedom 
[3]. The correct asymptotic behavior of this thermody- 
namical quantity can be reproduced using a careful and 
consistent implementation of the Pauli-Villars regular- 
ization both in NJL and the PNJL models [3] [4]. 

The extension to include the effect of the Polyakov 
loop can be done by introducing a background of static 
and diagonal (in the Polyakov gauge) temporal gluonic 
fields, A 4 = iA°, with A* 1 = 8^gA a ^- (where X" are the 
GellMann matrices) coupled to the quark fields by the 
covariant derivative: —> = + iA^. The quan- 



tities <p and <j) (which in the framework of the model 
are treated as classical field variables) are then given by 
the trace of the Polyakov loop operator and its adjoint 
denotes path ordering, N c is the number of colors 
and j3 is the inverse of the temperature): § = Ti c L/N c , 

= Tr c L f /N c , L = &>e fo ^4 . 

The transition between the confined-deconfined 
phases is driven by the temperature dependence of the 
additional pure gluonic term, the Polyakov potential, 
% ', for which several forms have been proposed. In [5] 
a polynomial form motivated by a Ginzburg-Landau 
ansatz was used and we will refer to it as . In [6] a 
form is used which we refer to as a i/ n , including a log- 
arithmic term inspired by the Haar measure of SU(N C ) 
group integration. In [7] an exponential term derived in 
the strong coupling expansion of the lattice QCD action 
is included, °lt 1 , and in [8] a form , combines the 
polynomial and logarithmic forms (for details see these 
references and our discussion of these in [4] where we 
also discuss the parametrization choice). 

Integrating the gap equations selfconsistently with the 
stationary phase equations, 

m u -M u = Gh u + ±h d h s + &h u hj + 
m d -M d = Gh d + ±hJi s + ^h d h 2 f +^h d , (1) 
m s - M s = Gh s + ±h u h d + S±h s hj + 

(where mf and Mf are the current and dynamical masses 
respectively), the thermodynamical potential is obtained 
[3] (T is the temperature, pi the chemical potential, G, 
K, gi and g2 are the coupling strengths of the of the 
NJL, 't Hooft, OZI-violating and non-violating 8-quark 



interactions): 



and for the medium contributions we obtain: 



£l(M f ,T,n,(j>,(j>) = 

=^ ( 4Ghj + Kh u h d h s + (h 2 f + 3g 2 hj 
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£ (j- 1 (Mj,T,n,<j>,<j))+C(T,n))+^ (0,0, r). 



f=u,d,s 



(2) 



The effect of the Polyakov loop was included straight- 
forwardly by noting that its phase enters the action as 
an imaginary chemical potential thus resulting in the fol- 
lowing generalizations of the usual expressions: 



1 



E p ,y,, T,0,0) = — Y,n q{ q)(E p ,n + i(A 4 ) u ,T), 



/+(") (E p ,n,T,<j>,<j>)=Logl\ 



-(E P TH)/T 



i=\ n q{q)( E P^ + l ( A ^ii^ T ) 

(3) 

Here E p = \~]?e\ 2 + M 2 and the (anti-)quark occupation 
numbers are give as usual by: n q ^ = (1 -f e ( £ pW)/ 7 ')- 1 



REGULARIZATION 

We compare the results obtained using the usual 3- 
dimensional cutoff in the momentum integrals with a 
Pauli-Villars covariant regularization with two subtrac- 
tions in the integrand. These are therefore specified by 
the action of the operators 



p 3D = ®(A-\rf\ 



(4) 



Pape = 1 - 1 



in the momentum integrand function. 

The vacuum and medium contributions can be sepa- 
rated and depending on the choice of the regularization 
procedure we obtain for the vacuum contributions, 

,vac(3D)/,.2\ 
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=A (2A 3 - VM 2 +A 2 (M 2 + 2A 2 ) J +M 4 ArcSinh— 



fac{PV) {M2) 



M 4 -A\ „ M 2 M 2 / 2 - M 2 \ 
= Ml + ^r)-— (A 2 +M 2 l nxr ), (6) 



(M 2 ,r, M ,0,0) 
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rf (w) (M 2 ,7>,0,0), 



Jo 



-pAp E 



n q M + n-qM _ n q o + nq p 
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Here we used the abbreviated notation with the sub- 
scripts M and denoting the quantities defined in (3) 
evaluated at those values for the mass (for the zero mass 
case we also set = = 1) [3] [4]. 

The constants of integration, C(T,jj,), resulting from 
the integration of the gap equations over the mass are 
chosen as to counterbalance the part stemming from the 
zero-mass limit of integration and are therefore given by: 

C 3 °(7» = - j\\pl\8\pl\ 2 T (/+ +/ ") , 
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C PV (T^)=- d\% 
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Note that while this may not be the standard way used 
to derive the 3D result the final result is equivalent to the 
usual one . 

A common variation of this procedure has been to re- 
move the cutoff of the medium parts as they are con- 
vergent. In the next section we compare results obtained 
using these four possibilities. 



RESULTS 

Here we only present some selected results which illus- 
trate the key points, for a more complete exposition and 
further details we refer to [4]. 

We tested the above mentioned four different regular- 
ization procedures with the four mentioned Polyakov po- 
tential forms, using parameter sets with and without eight 
quark interactions. The analysis of the normalized pres- 
sure difference at vanishing chemical potential which 
serves as a measure for the effective degrees of freedom, 
v(T) = (p(T)-p(0))/(7t 2 T 4 /90), reveals that the fail- 
ure to reach the Stefan-Boltzmann limit only happens in 
the case where we use the 3D cutoff everywhere (see Fig. 
la). 

The removal of the cutoff in the medium part results in 
the dynamical mass dropping below the current mass go- 
ing asymptotically to zero with the temperature increase 
(see Fig. lb). 

The removal of the cutoff leads to overshooting of the 
asymptotic solution dictated by the pure gluonic part in 
the case of potentials W 1 , W 11 and <% IV (in W 11 the 



logarithmic divergence prevents this from happening) as 
can be seen in Fig. 2a, 2b, 3a and 3b. This overshooting 
does not occur if the cutoff is kept (see Fig. 2b for the 
finite chemical potential case). 

The reason for these behaviours can be traced back 
to the asymptotic behaviour of the derivatives of /_i 
with respect to M, and <j> when T —> °°: the derivative 
with respect to the mass diverges upon removal of the 
cutoff whereas with cutoff it goes to zero; the derivative 
with respect to (0) diverges with T 4 upon the removal 
of the cutoff, the same order as the derivative of a tt ' , 
originating a deviation from the solution dictated by the 
latter. The divergence is lower with cutoff and in this 
case the asymptotic solution is dictated by the Polyakov 
potential. 




FIGURE 1. Temperature dependence ([r] = GeV) at \i = 
in the PNJL model without eight quark interactions and with 
Polyakov potential ty 1 of: the number of effective degrees of 
freedom, v(T), in la (black lines using PV and gray using 
3D, dashing denotes the removal of the cutoff in the medium 
contributions), dynamical mass of the quarks {[Mi] = GeV) 
using PV regularization ( thiner lines correspond to the removal 
of the cutoff in the medium). The results for the dynamical mass 
with 3D regularization are qualitatively similar. 




FIGURE 2. Temperature dependence of the Polyakov loop 
for the case withou 8q interactions and using a ^ 1 : a) at jl — 
(0 = (j>) the black lines refer to the use of PV and gray to 3D 
(dashed with only the vacuum regularized); b) <j> in full lines 
and <j> in dashed lines at = 100 MeV, black and gray lines 
refer to PV and 3D respectively (only the cases with cutoff 
everywhere are displayed). 




FIGURE 3. Same as in Figure 2a but using a ^/ lu without 8q 
interactions in the left panel and with 8q interactions in 
the right panel (same color code for the regularizations as in 
2a). 



CONCLUSIONS 

These qualitative features appear to be independent of 
the choice of parametrization (both of the quark interac- 
tions and the Polyakov potential) and are in fact a result 
of the regularization. The choice of Pauli- Villars regular- 
ization with the cutoff consistently kept over all contribu- 
tions achieves the best results for the studied quantities. 

ACKNOWLEDGMENTS 

This work has been supported in part by grants of 
Fundacao para a Ciencia e Tecnologia, FEDER, OE, 
SFRH/BPD/63070/2009, and Centra de Ffsica Computa- 
cional, unit 405. We acknowledge the support of the Eu- 
ropean Community-Research Infrastructure Integrating 
Activity Study of Strongly Interacting Matter (acronym 
HadronPhysics2, Grant Agreement No. 227431) under 
the Seventh Framework Programme of the EU. 



REFERENCES 

1. A. A. Osipov, B. Hiller, and J. da Providencia, Phys. Lett. 
B634, 48-54 (2006), hep-ph/0 50 80 5 8. 

2. A. A. Osipov, B. Hiller, A. H. Blin, and J. da Providencia, 
Annals Phys. 322, 2021-2054 (2007), hep-ph/ 
0607066' 

3. B. Hiller, J. Moreira, A. A. Osipov, and A. H. Blin, Phys. 
Rev. D81, 116005 (2010), 0812. 1532. 

4. J. Moreira, B. Hiller, A. A. Osipov, and A. H. Blin (2010), 
1008.0569. 

5. C. Ratti, M. A. Thaler, and W. Weise, Phys. Rev. D73, 
014019 (2006), hep-ph/0506234. 

6. S. Roessner, C. Ratti, and W. Weise, Phys. Rev. D75, 
034007 (2007), hep-ph/0 60 9281. 

7. K. Fukushima, J. Phys. G35, 104020 (2008), 
hep-ph/0806.0292. 

8. A. Bhattacharyya, P. Deb, S. K. Ghosh, and R. Ray (2010), 
1003.3337. 



